Role for Carboxylic Acid Moiety in NSAIDs: Favoring the Binding at Site II of Bovine Serum Albumin.

The molecular structures of nonsteroidal anti-inflammatory drugs (NSAIDs) vary, but most contain a carboxylic acid functional group (RCOOH). This functional group is known to be related to the mechanism of cyclooxygenase inhibition and also causes side effects, such as gastrointestinal bleeding. This study proposes a new role for RCOOH in NSAIDs: facilitating the interaction at the binding site II of serum albumins. We used bovine serum albumin (BSA) as a model to investigate the interactions with ligands at site II. Using dansyl-proline (DP) as a fluorescent site II marker, we demonstrated that only negatively charged NSAIDs such as ibuprofen (IBP), naproxen (NPX), diflunisal (DFS), and ketoprofen (KTP) can efficiently displace DP from the albumin binding site. We confirmed the importance of RCOO by neutralizing IBP and NPX through esterification, which reduced the displacement of DP. The competition was also monitored by stopped-flow experiments. While IBP and NPX displaced DP in less than 1 s, the ester derivatives were ineffective. We also observed a higher affinity of negatively charged NSAIDs using DFS as a probe and ultrafiltration experiments. Molecular docking simulations showed an essential salt bridge between the positively charged residues Arg409 and Lys413 with RCOO-, consistent with the experimental findings. We performed a ligand dissociation pathway and corresponding energy analysis by applying molecular dynamics. The dissociation of NPX showed a higher free energy barrier than its ester. Apart from BSA, we conducted some experimental studies with human serum albumin, and similar results were obtained, suggesting a general effect for other mammalian serum albumins. Our findings support that the RCOOH moiety affects not only the mechanism of action and side effects but also the pharmacokinetics of NSAIDs.

Insights into the interaction of human serum albumin with ionic liquids - Thermodynamic, spectroscopic and molecular modelling studies.

Human serum albumin (HSA) effectively binds different types of low-molecular-weight compounds and thus enables their distribution in living organisms. Recently, it has been reported that the protein-ligand interactions play a crucial role in bioaccumulation processes and provide an important sorption phase, especially for ionogenic compounds. Therefore, the binding interactions of such compounds with proteins are the subject of an ongoing interest in environmental and life sciences. In this paper, the influence of some counter-ions, namely [B(CN)4]- and [C(CN)3]- on the affinity of the [IM1-12]+ towards HSA has been investigated and discussed based on experimental methods (isothermal titration calorimetry and steady-state fluorescence spectroscopy) and molecular dynamics-based computational approaches. Furthermore, the thermal stability of the resulting HSA/ligand complexes was assessed using DSC and CD spectroscopy. As an outcome of the work, it has been ascertained that the protein is able to bind simultaneously the ligands under study but in different regions of HSA. Thus, the presence in the system of [IM1-12]+ does not disturb the binding of [C(CN)3]- and [B(CN)4]-. The presented results provide important information on the presence of globular proteins and some ionogenic compounds in the distribution and bioaccumulation of ILs in the environment and living organisms.

Cyclic Peptidic Furin Inhibitors Developed by Combinatorial Chemistry.

Furin is a human serine protease responsible for activating numerous physiologically relevant cell substrates and is also involved in the development of various pathological conditions, including inflammatory diseases, cancers, and viral and bacterial infections. Therefore, compounds with the ability to inhibit furin's proteolytic action are regarded as potential therapeutics. Here we took the combinatorial chemistry approach (library consisting of 2000 peptides) to obtain new, strong, and stable peptide furin inhibitors. The extensively studied trypsin inhibitor SFTI-1 was used as a leading structure. A selected monocylic inhibitor was further modified to finally yield five mono- or bicyclic furin inhibitors with values of K i in the subnanomolar range. Inhibitor 5 was the most active (K i = 0.21 nM) and significantly more proteolytically resistant than the reference furin inhibitor described in the literature. Moreover, it reduced furin-like activity in PANC-1 cell lysate. Detailed analysis of furin-inhibitor complexes using molecular dynamics simulations is also reported.

Repulsive Scaling Replica Exchange Molecular Dynamics in Modeling Protein-Glycosaminoglycan Complexes.

Glycosaminoglycans are long linear periodic anionic polysaccharides consisting of disaccharide units exhibiting different sulfation patterns forming a highly heterogeneous group of molecules. Due to their flexibility, length, high charge, and periodicity, they are challenging for computational approaches. Despite their biological significance in terms of the important role in various diseases (e.g., Alzheimer, cancer, SARS-CoV-2) and proper cell functioning (e.g., proliferation, maturation), there is a lack of effective molecular docking tools designed specifically for glycosaminoglycans due to their challenging physical-chemical nature. In this chapter we present protocols for the Repulsive Scaling Replica Exchange Molecular Dynamics (RS-REMD) methods to dock glycosaminoglycans with both implicit and explicit solvent models implemented. This novel molecular dynamics-based replica exchange technique should help to elevate our current knowledge on the complexes and interactions between glycosaminoglycans and their protein receptors.

Binding of heparan sulfate to human cystatin C modulates inhibition of cathepsin L: Putative consequences in mucopolysaccharidosis.

Mucopolysaccharidoses (MPS) are a group of rare lysosomal storage diseases characterized by glycosaminoglycan (GAG) accumulation causing progressive multi-organs dysfunction and ultimately severe cardio-respiratory damages. Human cystatin C (hCC), a potent inhibitor of cysteine cathepsins, plays an important role in respiratory diseases. However, its regulation remained unknown in MPS. Herein, elevated hCC levels were measured in respiratory specimens from MPS-I, -II, and -III patients and were significantly correlated with severe respiratory symptoms (rs = 0.7173). Heparan sulfate (HS), a prominent GAG, dampened its inhibitory activity toward cathepsin L in a dose-dependent manner. HS and HS-oligosaccharides bound tightly hCC, in combination with a secondary structure rearrangement. Molecular modeling studies identified three HS binding regions in hCC, including the N-terminus, which is crucial in the inhibition of cathepsins. Impairment of inhibitory potential of hCC may reflect abnormal regulation of proteolytic activity of cathepsin L in lung, ultimately contributing to the severity of MPS.

Molecular dynamics-based descriptors of 3-O-Sulfated Heparan sulfate as contributors of protein binding specificity.

Glycosaminoglycans are linear periodic and anionic polysaccharides found in the extracellular matrix, involved in a range of key biochemical processes as a result of their interactions with a variety of protein partners. Due to the template-less synthesis, high flexibility and charge of GAGs, as well as the multipose binding of GAG ligands to receptors, the specificity of GAG-protein interactions can be difficult to elucidate. In this study we propose a set of MD-based descriptors of unbound Heparan Sulfate hexasaccharides that can be used to characterize GAGs and explain their binding affinity to a set of protein receptors. With the help of experimental data on GAG-protein binding affinity, we were able to further characterize the nature of this interaction in addition to providing a basis for predictor functions of GAG-protein binding specificity.

Molecular Dynamics Approaches Dissect Molecular Mechanisms Underlying Methylene Blue-Glycosaminoglycan Interactions.

Glycosaminoglycans (GAGs) are a class of periodic anionic linear polysaccharides involved in a number of biologically relevant processes in the extracellular matrix via interactions with various types of molecules including proteins, peptides and small organic molecules. The metachromatic dye methylene blue (MB) is a GAG binding agent. This molecule possesses a tricyclic, monocationic phenothiazine ring system, while the terminal methyl groups attached to the nitrogen atoms bear the most positive charges of the cation and, therefore, represent potential binding sites for negatively charged GAGs. In this study, we rigorously explored molecular mechanisms underlying these interactions for several GAG types: heparin, heparan and chondroitin sulfates. We found that GAG-MB interactions are predominantly electrostatically driven, with the particularly important role of sulfate groups. MB oligomeric stack formation was favored in the presence of GAGs. Furthermore, the impact of MB binding on the conformation of GAGs was also evaluated. The novel results allow for better quantitative analytics of GAG composition in the studied biochemical systems using MB dye as a GAG-specific marker. Our data add to the knowledge on small molecule-GAG interactions and could be potentially useful for novel developments in drug design and putative disease therapies in which GAGs are involved.

Induced circular dichroism as a tool to monitor the displacement of ligands between albumins.

The induction of chirality in a ligand can be a powerful analytical tool for studying protein-ligand interactions. Here, we advanced by applying the technique to monitor the inversion of the induced circular dichroism (ICD) spectrum when ligands move between human and bovine serum albumin proteins (HSA and BSA). ICD experiments were performed using dimers of methyl vanillate (DVT) and vanillin (DVN). The sign and spectra shape were dependent on the albumin type. DVN presented a positive maximum in 312 nm when complexed with HSA and a negative one in BSA. It was possible to induce and follow the time-dependent displacement of the ligand from BSA (2.2 × 106 M-1) to HSA (6.6 × 105 M-1) via ICD inversion. The Molecular Mechanics Generalized Born Surface Area approach was used to calculate the binding free energy of the conformers, and a dissociation pathway for each system was proposed using Umbrella Sampling calculations. Four energy minima dihedral angle conformers were identified, and the corresponding CD spectra were calculated using the quantum chemistry approach. Then, weighted spectra for the conformationally accessible conformers were obtained based on each conformer's Boltzmann probability distribution. In conclusion, the methodology described in the manuscript might be helpful in monitoring the movement of ligands between proteins that they bind.

Advanced Molecular Dynamics Approaches to Model a Tertiary Complex APRIL/TACI with Long Glycosaminoglycans.

Glycosaminoglycans (GAGs) are linear anionic periodic polysaccharides participating in a number of biologically relevant processes in the extracellular matrix via interactions with their protein targets. Due to their periodicity, conformational flexibility, pseudo-symmetry of the sulfation pattern, and the key role of electrostatics, these molecules are challenging for both experimental and theoretical approaches. In particular, conventional molecular docking applied for GAGs longer than 10-mer experiences severe difficulties. In this work, for the first time, 24- and 48-meric GAGs were docked using all-atomic repulsive-scaling Hamiltonian replica exchange molecular dynamics (RS-REMD), a novel methodology based on replicas with van der Waals radii of interacting molecules being scaled. This approach performed well for proteins complexed with oligomeric GAGs and is independent of their length, which distinguishes it from other molecular docking approaches. We built a model of long GAGs in complex with a proliferation-inducing ligand (APRIL) prebound to its receptors, the B cell maturation antigen and the transmembrane activator and calcium modulator and cyclophilin ligand interactor (TACI). Furthermore, the prediction power of the RS-REMD for this tertiary complex was evaluated. We conclude that the TACI-GAG interaction could be potentially amplified by TACI's binding to APRIL. RS-REMD outperformed Autodock3, the docking program previously proven the best for short GAGs.

Prediction of protein assemblies, the next frontier: The CASP14-CAPRI experiment.

We present the results for CAPRI Round 50, the fourth joint CASP-CAPRI protein assembly prediction challenge. The Round comprised a total of twelve targets, including six dimers, three trimers, and three higher-order oligomers. Four of these were easy targets, for which good structural templates were available either for the full assembly, or for the main interfaces (of the higher-order oligomers). Eight were difficult targets for which only distantly related templates were found for the individual subunits. Twenty-five CAPRI groups including eight automatic servers submitted ~1250 models per target. Twenty groups including six servers participated in the CAPRI scoring challenge submitted ~190 models per target. The accuracy of the predicted models was evaluated using the classical CAPRI criteria. The prediction performance was measured by a weighted scoring scheme that takes into account the number of models of acceptable quality or higher submitted by each group as part of their five top-ranking models. Compared to the previous CASP-CAPRI challenge, top performing groups submitted such models for a larger fraction (70-75%) of the targets in this Round, but fewer of these models were of high accuracy. Scorer groups achieved stronger performance with more groups submitting correct models for 70-80% of the targets or achieving high accuracy predictions. Servers performed less well in general, except for the MDOCKPP and LZERD servers, who performed on par with human groups. In addition to these results, major advances in methodology are discussed, providing an informative overview of where the prediction of protein assemblies currently stands.

Modeling protein structures with the coarse-grained UNRES force field in the CASP14 experiment.

The UNited RESidue (UNRES) force field was tested in the 14th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP14), in which larger oligomeric and multimeric targets were present compared to previous editions. Three prediction modes were tested (i) ab initio (the UNRES group), (ii) contact-assisted (the UNRES-contact group), and (iii) template-assisted (the UNRES-template group). For most of the targets, the contact restraints were derived from the server models top-ranked by the DeepQA method, while the DNCON2 method was used for 11 targets. Our consensus-fragment procedure was used to run template-assisted predictions. Each group also processed the Nuclear Magnetic Resonance (NMR)- and Small Angle X-Ray Scattering (SAXS)-data assisted targets. The average Global Distance Test Total Score (GDT_TS) of the 'Model 1' predictions were 29.17, 39.32, and 56.37 for the UNRES, UNRES-contact, and UNRES-template predictions, respectively, increasing by 0.53, 2.24, and 3.76, respectively, compared to CASP13. It was also found that the GDT_TS of the UNRES models obtained in ab initio mode and in the contact-assisted mode decreases with the square root of chain length, while the exponent in this relationship is 0.20 for the UNRES-template group models and 0.11 for the best performing AlphaFold2 models, which suggests that incorporation of database information, which stems from protein evolution, brings in long-range correlations, thus enabling the correction of force-field inaccuracies.

The potential role of glycosaminoglycans in serum amyloid A fibril formation by in silico approaches.

Serum amyloid A (SAA) is actively involved in such pathological processes as atherosclerosis, rheumatoid arthritis, cancer and Alzheimer's disease by its aggregation. One of the factors that can attenuate its aggregation and so affects its physiological role is its interactions with glycosminoglycans (GAGs), linear anionic periodic polysaccharides. These molecules located in the extracellular matrix of the cell are highly variable in their chemical composition and sulfation patterns. Despite the available experimental evidence of SAA-GAG interactions, no mechanistic details at atomic level have been reported for these systems so far. In our work we aimed to apply diverse computational tools to characterize SAA-GAG complexes formation and to answer questions about their potential specificity, energetic patterns, particular SAA residues involved in these interactions, favourable oligomeric state of the protein and the potential influence of GAGs on SAA aggregation. Molecular docking, conventional and replica exchange molecular dynamics approaches were applied to corroborate the experimental knowledge and to propose the corresponding molecular models. SAA-GAG complex formation was found to be electrostatics-driven and rather unspecific of a GAG sulfation pattern, more favorable for the dimer than for the monomer when binding to a short GAG oligosaccharide through its N-terminal helix, potentially contributing to the unfolding of this helix, which could lead to the promotion of the protein aggregation. The data obtained add to the specific knowledge on SAA-GAG systems and deepen the general understanding of protein-GAG interactions that is of a considerable value for the development of GAG-based approaches in a broad theurapeutic context.

Computational insights into heparin-small molecule interactions: Evaluation of the balance between stacking and non-stacking binding modes.

Glycosaminoglycans (GAGs), anionic periodic linear polysaccharides, are involved in a manifold of key biochemical processes ongoing in the extracellular matrix via establishing direct intermolecular interactions with diverse classes of biopolymers as well as with bioactive small molecules. Due to their acidic nature, they are capable of binding positively charged ligands, which, in turn could affect their binding with protein and peptide targets, modulating a number of physiologically important signaling pathways. Therefore, it is of great significance to improve our understanding on the molecular basis underlying GAG-small molecule interactions. In this study, we applied in silico approaches (molecular dynamics and free energy calculations) complemented with circular dichroism and absorption spectroscopy to characterize the complex formation between heparin, one of the principal members of GAG family, and twenty different cationic ligands including therapeutic drugs, alkaloids and organic dyes. In particular, the oligomerization propensity of ligands prior to heparin binding, binding free energy parameters, effects of the ionic strength are rigorously described. Based on the performed analysis, the ligands are classified into three main groups depending on their heparin binding and oligomerization properties. The computational data agree and provide rationale for the corresponding experimental findings, contributing to the general knowledge of the physico-chemical nature of ligand-GAG intermolecular interactions.

Evaluation of replica exchange with repulsive scaling approach for docking glycosaminoglycans.

Glycosaminoglycans (GAGs), long linear periodic anionic polysaccharides, are key molecules in the extracellular matrix (ECM). Therefore, deciphering their role in the biologically relevant context is important for fundamental understanding of the processes ongoing in ECM and for establishing new strategies in the regenerative medicine. Although GAGs represent a number of computational challenges, molecular docking is a powerful tool for analysis of their interactions. Despite the recent development of GAG-specific docking approaches, there is plenty of room for improvement. Here, replica exchange molecular dynamics with repulsive scaling (REMD-RS) recently proved to be a successful approach for protein-protein complexes, was applied to dock GAGs. In this method, effective pairwise radii are increased in different Hamiltonian replicas. REMD-RS is shown to be an attractive alternative to classical docking approaches for GAGs. This work contributes to setting up of GAG-specific computational protocols and provides new insights into the nature of these biological systems.

Acute phase α1-acid glycoprotein as a siderophore-capturing component of the human plasma: A molecular modeling study.

Siderophores are ferric ion-specific organic compounds that are used by bacteria and fungi to secure their iron supply when infecting target organisms. There are a few proteins in the human body, named siderocalins, which bind these important virulence factors and so starve microorganisms of iron. In this study, we analyzed in silico if serum α1-acid glycoprotein (AAG), the major acute phase lipocalin component of the human plasma, could functionally belong to this group. The real biological function of AAG is elusive and its concentration substantially increases in response to pathological stimuli, including bacterial infections. We computationally evaluated the potential binding of nine microbial siderophores into the ß-barrel cavity of AAG and compared the results with the corresponding experimental data reported for siderophore-neutrophil gelatinase-associated lipocalin complexes. According to the results, petrobactin and Fe-BisHaCam are putative candidates to be recognized by this protein. It is proposed that AAG may function as a siderophore capturing component of the innate immune system being able to neutralize bacterial iron chelators not recognized by other siderocalins.

Computational insights into the role of calcium ions in protein-glycosaminoglycan systems.

Glycosaminoglycans (GAGs) are anionic, periodic, linear polysaccharides which are composed of periodic disaccharide units. They play a vital role in many biological processes ongoing in the extracellular matrix. In terms of computational approaches, GAGs are very challenging molecules due to their high flexibility, periodicity, predominantly electrostatic-driven nature of interactions with their protein counterparts and potential multipose binding. Furthermore, the molecular mechanisms underlying GAG-mediated interactions are not fully known yet, and experimental techniques alone are not always sufficient to gain insights into them. The aim of this study was to characterize protein-ion-GAG complexes for the systems where ions are directly involved in GAG binding. Molecular docking, molecular dynamics and free energy calculation approaches were applied to model and rigorously analyse the interactions between annexins (II and V), calcium ions (Ca2+) and heparin (HP). The computational data were examined and discussed in the context of the structural data previously reported by the crystallographic studies. The computational results confirm that the presence of Ca2+ has a tremendous impact on the annexin-HP binding site. This study provides a general computational pipeline to discover the complexity of protein-GAG interactions and helps to understand the role of ions involved at the atomic level. The limitations of the applied protocols are described and discussed pointing at the challenges persisting in the state-of-the-art in silico tools to study protein-ion-GAG systems.

Mind Your Dye: The Amyloid Sensor Thioflavin T Interacts with Sulfated Glycosaminoglycans Used To Induce Cross-ß-Sheet Motifs.

Thioflavin T (ThT) is a commonly employed fluorescence probe for sensing amyloid fibrils. These highly ordered, insoluble nanostructures colocalize with sulfated glycosaminoglycans (GAGs) being abundant in the extracellular matrix and on the outer surface of cell membranes. To elucidate the positive impact of GAGs on amyloidogenesis, they are frequently used to promote fibril formation in vitro, which is detected by the enhanced fluorescence of ThT. The polyanionic nature and the affinity of GAGs to basic compounds predict that cationic ThT molecules may also bind to them, in addition to cross-ß-structures formed in the reaction medium. By means of circular dichroism (CD) and absorption spectroscopy, this study examined the heparin and chondroitin sulfate binding of ThT. The large blue shift of the absorption peak indicated a card-pack type oligomerization of the dye molecules along the linear GAG chains. The strong exciton couplet observed in the CD spectra implies the left-handed, helical arrangement of GAG-associated oligomers of the dye. The decisive contribution of ionic forces for the binding was illustrated by sodium-ion-provoked dissociation of dye-GAG complexes. In silico analysis was performed to complement the experimental findings and to contribute to the understanding of potential molecular mechanisms underlying ThT-GAG interactions. ThT can be considered as an inert component in GAG-induced amyloid assays but only if the experiments are correctly designed.

Disulfide-Linked Peptides for Blocking BTLA/HVEM Binding.

Immune checkpoints are crucial in the maintenance of antitumor immune responses. The activation or blockade of immune checkpoints is dependent on the interactions between receptors and ligands; such interactions can provide inhibitory or stimulatory signals, including the enhancement or suppression of T-cell proliferation, differentiation, and/or cytokine secretion. B-and T-lymphocyte attenuator (BTLA) is a lymphoid-specific cell surface receptor which is present on T-cells and interacts with herpes virus entry mediator (HVEM), which is present on tumor cells. The binding of HVEM to BTLA triggers an inhibitory signal which attenuates the immune response. This feature is interesting for studying the molecular interactions between HVEM and BTLA, as they may be targeted for novel immunotherapies. This work was based on the crystal structure of the BTLA/HVEM complex showing that BTLA binds the N-terminal cysteine-rich domain of HVEM. We investigated the amino acid sequence of HVEM and used molecular modeling methods to develop inhibitors of the BTLA/HVEM interaction. We synthesized novel compounds and determined their ability to interact with the BTLA protein and inhibit the formation of the BTLA/HVEM complex. Our results suggest that the HVEM (14-39) peptide is a potent inhibitor of the formation of the BTLA/HVEM protein complex.

NMR and crystallographic structural studies of the extremely stable monomeric variant of human cystatin C with single amino acid substitution.

Human cystatin C (hCC), a member of the superfamily of papain-like cysteine protease inhibitors, is the most widespread cystatin in human body fluids. This small protein, in addition to its physiological function, is involved in various diseases, including cerebral amyloid angiopathy, cerebral hemorrhage, stroke, and dementia. Physiologically active hCC is a monomer. However, all structural studies based on crystallization led to the dimeric structure formed as a result of a three-dimensional exchange of the protein domains (3D domain swapping). The monomeric structure was obtained only for hCC variant V57N and for the protein stabilized by an additional disulfide bridge. With this study, we extend the number of models of monomeric hCC by an additional hCC variant with a single amino acid substitution in the flexible loop L1. The V57G variant was chosen for the X-ray and NMR structural analysis due to its exceptional conformational stability in solution. In this work, we show for the first time the structural and dynamics studies of human cystatin C variant in solution. We were also able to compare these data with the crystal structure of the hCC V57G and with other cystatins. The overall cystatin fold is retained in the solute form. Additionally, structural information concerning the N terminus was obtained during our studies and presented for the first time. DATABASE: Crystallographic structure: structural data are available in PDB databases under the accession number 6ROA. NMR structure: structural data are available in PDB and BMRB databases under the accession numbers 6RPV and 34399, respectively.

Characteristics of C-terminal, ß-amyloid peptide binding fragment of neuroprotective protease inhibitor, cystatin C.

Cystatin C originally identified as a cysteine proteases inhibitor has a broad spectrum of biological roles ranging from inhibition of extracellular cysteine protease activities, bone resorption, and modulation of inflammatory responses to stimulation of fibroblasts proliferation. There is an increasing number of evidence to suggest that human cystatin C (hCC) might play a protective role in the pathophysiology of sporadic Alzheimer's disease. In vivo and in vitro results well documented the association of hCC with Aß and the hCC-induced inhibition of Aß fibril formation. In our earlier work, using a combination of selective proteolytic methods and MS spectroscopy, C-terminal fragment hCC(101-117) was identified as the Aß-binding region. The fragment of Aß peptide responsible for the complex formation with hCC was found in the middle, highly hydrophobic part, Aß(17-24). Structures and affinities of both Aß and hCC binding sites were characterized by the enzyme-linked immunosorbent assay-like assay, by surface plasmon resonance, and by nano-ESI-FTICR MS of the hCC-Aß-binding peptide complexes. In the in vitro inhibition studies, the binding cystatin sequence, hCC(101-117), revealed the highest relative inhibitory effect toward Aß-fibril formation. Herein, we present further studies on molecular details of the hCC-Aß complex. With Ala substitution, affinity experiments, and enzyme-linked immunosorbent assay-like assays for the Aß-binding fragment, hCC(101-117), and its variants, the importance of individual amino acid residues for the protein interaction was evaluated. The results were analyzed using hCC(101-117) nuclear magnetic resonance structural data with molecular dynamics calculations and molecular modeling of the complexes. The results point to conformational requirements and special importance of some amino acid residues for the protein interaction. The obtained results might be helpful for the design of low molecular compounds modulating the biological role of both proteins. Copyright © 2016 John Wiley & Sons, Ltd.